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CHOICE OF PROFILE FOR THE WINGS OF AN AIRPLANE. * 
By A. Toussaint and E. Carafoli. 

Simplified General Method for Drawing Airplane Profiles 

The above method was useful for obtaining a better under- 
standing of the problem,* It leads, nevertheless, to a somewhat 
laborious drawing which becomes quite complicated when we take 
a transformation- function having terms of a high degree. 

Besides , it is extremely difficult to determine by that 
method the parameters which govern the profile form, with the 
object of obtaining "families of profiles," that is^ profiles 
the form of which evolves in a continuous manner by modification 
of one of the characteristics, such as 0 rao , the maximum rel- 
ative thickness; e/i , the suitable distribution of the relative 
thickness along the chord, etc. 



*From L 1 Aeronaut ique, # January, 1928. For Part I see N.A.C..A. 
Technical Memorandum *No# 468 • 



N.A.C.A. Technical Memorandum .No. 469 



2 



The method given below overcomes these difficulties. It 
depends on the following considerations (See Comptes. Rendus de 
I'Academie des Sciences, "Notes" de M. Carafoli, October 24,. 
1927, and November 14 and 28, 1927. 

We have seen that, by changing the axes, the transforma- 
tion functioir can always be put into the form 

z= ^f + S ^ (12') 

in which c 2 .is a real quantity. 

Let Ox and Oy be a system of axes, M the. center of the 
generating circle passing through B x , and B* the point of 
zero velocity which should correspond to the profile tip' (Fig. 
11). It can be shown that the auxiliary circle VL 9 corre- 
sponding to the partial transformation 

z. = 2i or z t = c 2 

passes through the point C 1 so that 0B 1 x 00 x = cs, and 
the center M l is located at the intersection of the line 0M X 
(symmetrical to OM with respect to Oy) with the line C x M x 
parallel to B X M (Fig. 8). 

The position of B* on the generating circle is defined by 
the value of the moment coefficient C mo , according to the for- 
mula 

>o = 8TT fa T ( 10! ) 
in which t, the angle of B ! M with Ox, corresponds to 



N.A/C.A. Technical Memorandum No. 469 3 

.(0 - 7) of the preceding- notations (if* MB*0 = 0 and B l OB 1 
= 7) and I is the profile chord, the value of which is approxi- 
mately 2 (0B X + 00 l ). 

It follows from equation (10* ), t being expressed in de- 
grees, that 

„ 0m o 180 4(08,4-00, f 

8tt tT~ OB, x OCi " ° m o 36,4 " ^7 - 

OB, 

= 36.4 C mQ 

i 7 

For the point B 1 U = - 0B f e~ ), which must correspond 

to the trailing edge of the profile, we should have 

n 

S - — - 2 n x n 

(tf) B .° s c r " -o. (3 5 ) 

If the point C ! corresponding to B ! for the transforma- 
tion z x = be on the auxiliary circle 1^ , then (£ - ~\ 
represents CtB'e 10 (Fig. 11), from which we derive 

^B'e iG -(J 2 = 0 

or * (36) 
y 2 - nipt© 1 

(-l) n '(0B.) n e- lrt¥ 

If, for example, the transformation function be 



= » ; = si + ffB + fa + 

5 \(P f 1 ^ 1 



the above formula can then be written 
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(- 1)P (0B') P e- 1 ^ (- l) q (0B') q e- 1 ^ 

(36') 



rXr = C'B'e ia 



(.1)* <OB') r e l * 

so that the three terms of the first member constitute a polygon 
O^IHB 1 the geometrical resultant of which is B'C (Fig. 12). 

Conversely, each polygons CIHB' would correspond to three 
terms p, q, r of the transformation function, which are deter- 
mined "by means of the formula (36* ) 

(- 1) P (0B') P 

and so on. ' 

In particular, the polygon could be reduced to any number 
of sides aligned on C'B 1 (as, for example, CI'H'B 1 ). 
In practice, it even suffices- to consider a single supplementary 
term of any degree n, in which case we have 

(- l) n (0B') n 



♦We have already pointed out the advantage of the case in which 
the only complementary term is of the degree n = 3, or 



C2 /#,\ 
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We should thus add to the two terms, which give us the point 
P, the segment P P n of the modulus 

c'B 1 OB 1 n 



Tn 



= (- D J 



n OP 1 



and of amplitude 



qp n = a - n7 - n9 



In order to compute the velocity W p at the point' P n of 
the profile, it suffices to construct the vector 

m 

n 



Q'R 1 = 



n x 



7" 



n 



= (- i) n CB- (2|i) 



starting from the corresponding point Q f with direction cp n 
and magnitude P'R 1 , and we get, as before, 



W P - 27 "MB" X P¥ 



(25 ! ) 



Utilizing these principles, the profile design would involve 
the following operations, which are entirely geometrical ajid re- 
quire no special knowledge of the problem. 

The simplified method of drawing a profile with a pointed 
trailing edge includes the following operations: 

1. Draw the two axes Ox and Oy and lay out the length 
0B X along the negative direction of Ox. OB^ is generally 
taken as one-fourth, the chord of the profile to be obtained 
(Fig. 11). 

2 V . Choose OCi or B 1 C 1 determining the point C x , the 
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position of which relative to B 1 characterizes the thickness 
distribution along the profile chord* In the case represented by 
the figure, we have B 1 C 1 positive or negative according to 
whether C x is on the right or on the : left of B x .. When B x C x 
is positive, the profile is less sharp at the rear and conversely. 
In practice, we take 

t 0.05. 

0 B x 

3c Choose 0M o which characterizes the apparent mean cam- 
ber of the profile (by analogy with the case of Joukowski pro- 
files). For profiles with small C mo it is convenient to take 

™z < 0 . 10 . 

' 4. Choose the center M of the generating circle on the 
straight line B 1 M 0 in such a manner that M 0 M would satisfy 
the following formula in terms of the maximum relative thickness 
e/l we wish to obtain 

M 0 M ± 0,4 B 1 C 1 0*77 1 



0B 1 1 - 0.77 f 



The minus sign is taken when B 1 C x is positive and vice versa. 
5. Draw the generating circle with the center M and the 



radius MB X . 



6. Draw the auxiliary circle corresponding to the partial 



0 
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transformation* 

z 1 = or z£ = c 2 = 0B X x 0C X . 

This circle will then pass through C x . Its center 1^ 
will be at the intersection: of the line 0M X (symmetrical to 
x Oil with respect to Oy) with the line parallel to B X M. ' 

?• Carry out the operations corresponding to the first two 
terms of the function, by geometric summation of the vectors 
OP 1 and OQ* of equal amplitudes and opposite signs (Trefftz, 
geometrical construction). 

When G 1 is inside the generating circle, the points P 
define a profile with rounded trailing edge, but, when C x 
is outside of the circle, the contour of the points P forms 
a figure 8 with a small bulge toward the rear. 

8. Define, on the generating circle, the point of zero ve- 
locity B< , with the object of obtaining a 0^ fixed a priori . 
For the latter, we draw E'M making an angle r with Ox, such 
that t° = 36.4 C mo . The point C 1 , on the auxiliary circle 

M x , corresponds to the point B* , so that 

B 1 0B X = ^OC* = y 

9. Construct the representative segments of the other terms 
of the transformation function. 

Since, in practice, it is important to take a single com- 
plementary term of the degree n = 3, the segment to be con- 



N.A.C.A. Technical Memorandum No. 469 8 

struct ed, starting from each of the points P, will have the mag- 
nitude 



£k C'B 1 /OB'N 3 

£3 ~ " 3 VPP'/ 



and the direction 

<p = a _ 37 - 39 

Remarks . - The above method defines, without ambiguity, the 
parameters which govern the evolution of the profile form* These 
parameters are 

. DM*. MpM». 0.4 B^ . \G 1 . 
m o ' OBi ' 0B 1 ' OBi ' 

and the degree n of the complementary term. 

The knowledge of these parameters permits the definition of 

» 

"profile groups" of evolutive forms and characteristics. For ex- 
ample: 

(a) Groups with variable C mQ , with constant apparent 
camber and relative thickness and otherwise general form charac- 
terized by the same B 1 Q 1 /Q% 1 and the same degree n. 

(b) Groups with variable apparent camber, with constant 
C mo and relative thickness, as well as same B^/OBx and n. 

(c) Groups with variable thickness distribution (by varia- 
tion of B 1 G x /0B 1 ) 9 all other characteristics remaining constant, 

(d) Groups with variable relative thickness, with all other 
parameters constant. 
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(e) Groups with variable n, and all other characteristics 
constant. In practice, there seems to "be little interest, in 
the present state of the problem, in varying n and n = 3 is 
therefore adopted. 

Remarks .- The preceding general method also permits the 
dr a wing of profiles with rounded trailing edges (Comptes Rendus, 
Vol. 185, p. 842, Oct. 24, 1927). It can be adapted to the partic- 
ular case of complementary terms leading to corrective terms of 
constant modulus ( Comptes Rendus, Vol. 185, p. 1014, November, 1927). 

By extension we can pass from a pointed profile to a profile 
.. with a reflexed trailing edge, (Comptes Rendus, Vol. 185, p. 1189, 
Nov. 28, 1937. 

Aerodynamic Characteristics of Theoretical Profiles 

The aerodynamic characteristics of theoretical profiles re- 
sult from the application of the same conformal transformation 
to the sustaining flow around the generating circle, according 
to Joukowski's hypothesis. 

It is known that this hypothesis consists in determining 
the circulation V such that the point of the profile in the z 
pl'ane corresponds to a point of zero velocity on the generating 
circle in the £ plane. 

Under these conditions, the aerodynamic characteristics of 
the theoretical profile, corresponding to the case of a wing of 
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infinite span, will be given by the following formulas: 
Lift coefficient 

0 Z Z 8tt Jt (a ) (17) 
Z 

Moment coefficient at the angle of zero lift 
in which 

a is the angle of incidence with respect to OB 1 ; 

(3 is the amplitude of the first axis with respect to OB 1 ; . 

7 is the amplitude of the second axis with respect to OB 1 ; 

c 2 is the modulus of the parameter x x ; 

a is the radius of the generating circle;' 

I is the profile chord, measured from the drawing, after 
it is completed. (It is found, moreover, that the value thus 
obtained differs but little from the approached value.) 

In order to get the moment coefficient with reference to 
the leading edge of the profile, it is sufficient to determine 
the position: of the profile focus F. The latter is found on a 
line MF making the angle 2 (Y - 0 ) with the first axis ( MB f ) • 

The distance MF can then be calculated by the formula 

MF = — (19) . 

a 

Knowing thus the position of the focus F, we then measure ■ 
the distance FA as fax as the leading edge and obtain: 
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Cm 



(A) 



= °ra (F) ▼ 



or 



(30) 



° m (A) = wra 



= C m . 4- C z M 



In general M. is approximately 6,25. 



Diagrams of Aerodynamic Pressures on a Theoretical Profile 

The velocity at a given point on the profile being W p> 
the pressure coefficient can be defined by the expression 



Op = 1 



(21 } 



V being the aerodynamic velocity. 

The velocity W p is calculated by the general formula 



Wp = W c 



d£ 



(22) 



in which W c is the velocity corresponding to the point on the 
generating circle* 

For all points P» of the generating circle (Fig* 13), we 

have 

_ Dint mm 

. (23) 



W c = 2V £1^- = 8V » 



MP 

P ! D f being the length of the perpendicular dropped from P 1 to 
the line V which passes through B 1 • 
Moreover, we have 



9 rfi « 
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V 



1 


dz 









t - r - 



2x 2 eT^f ruc n e-U+IW 



r 3 C 

£' = OP' • 

Considering the point Q' conjugate of P l on the auxil- 
iary circle M x (Fig. 14), we have 



12 



(24) 



or 



If we introduce 



2x, e 



-3iY 



V 



of the amplitude cp 2 previously calculated, we get 



P'R'j. = P'Q« - Q'R'j = 

Similarly 



c s 2x 2 e 



-ai*Y 



3x 3 e 



-4lY 



r 



t 3 



with the amplitude cp 3 or (cp 3 + V) and we get 

2* a e-si'Y j 3^ e"^ 



P'R' S = 



V 



and so on, as far as the point R' (n- i) corresponding to the 
term in oc n . 

. It must "be noted that Q'R^ = 2 times the positional cor- 
rection previously calculated for the term in oc 3 . In like man- 
ner R» ( n _ 1 )R I ( n _ z y=H times' the correction calculated for the 
term in oc n . 

The construction of the segment P'R« ( n-l ) utilizes there- 
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fore the elements of the profile design. 

i- 

Finally, for each point P, consequent on P f , we have 

% - * ^ ifrsro • < 25 > 

from which C p is derived. 

The factors are calculable, once for all, for all 

p,R, U-i) 
points of the profile. 

The factors PD 1 vary with the angle of attack a for each 
of the points P ! . 

Aerodynamic Characteristics of a Wing of Finite Span, 
with a Theoretical Profile 

1. Case of rectangular wings 

The aerodynamic characteristics of a wing having a finite 
span are deduced from the aerodynamic characteristics of the pro- 
file of a wing of infinite span, by allowing for the modifications 
brought about in the theoretical flow by limiting the span. 

Prandtl's theory teaches us that these modifications are in- 
terpreted £n practice by the fact that: "The effective angle of 
attack at all sections of the finite wing is equal to the real 
angle of attack diminished by the. induced angle of the section: in 
question." This induced angle depends essentially on the law of 
lift distribution along the spaa. 

For rectangular monoplane wings the simplifying hypothesis 
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of an elliptical distribution of the lift is frequently used. 

a 

It can be considered as/suf fdciently practical approximation when 
we wish to calculate the mean aerodynamic .characteristics for 
the whole wing. If we wish to analyze what take's place in a. 
particular wing section, it is convenient to assume a lift dis~ 
tribution along the wing. span. In this connection, we have pro- 
posed to adopt the simplified approximate solution of Fuchs*, 
because the theoretical solution! of Betz necessitates very labor- 
ious calculations. The former consists in assuming that the law 
of distribution of the circulation! or of the lift on a rectangu- 
lar wing cf span A = £ ' ' is' 'represented by the formula 

1 + A /2oc\ 3 ' 
2A + 3tt VTW ' 

Figure 12 represents the comparison of the lift distribu- 
tions for a rectangular wing of aspect ratio 5 and £or the same 
mean lift of 1. 

The curve EEE corresponds to an elliptical distribution. 
The curve BBB corresponds to the evolutive distribution, as ex- 
pressed by formula (26). 

The straight line KR corresponds to a rectangular distri- 
bution:. 

Lastly, the dotted line represents the distribution- experi- 
mentally determined for a mean angle of attack. 

A sufficient agreement is noted between the theoretical 

distribution obtained by formula (26) and the experimental dis- 
*A. Tous saint and E. Carafoli, "Verifications experiment ales de 

nooJ he ° rie deS ailes sust sntatrices, « L'Aerophile, April 1-15, 
1927* 



» 
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tribution, except at the wing tips, where the latter exhibits 

abnormal deviations due to disturbances of the lateral edges. 

In adopting this distribution, approximate but at the same 

time developed with A, one easily finds that the mean lift 

coefficient C z on the wing is given by the formula 

a i 

1 + TT 1 



>Z U) z (°0 * 1 + 3n# (1 - 



J 



in which 



and 



= 0 , I -^4 (27) 

0 = stt^ (a + 0) (17) 



'•(co) i 
A 



a, = 



2A + Sir 



As before, a is the radius of the generating circle. It is 
sometimes expressed by an approximate formula in terms of the 
profile chord I. Since we are dealing with theoretical profiles, 
it is not necessary to resort to this approximation. 
The expression for Cz( A ) ma y also be written 



(27«)" 



°«(A) " A (A) (<X + P) 

with a , 

l + x*- 

A = 8 - tt a 5 ^ _ 

(A) I ^ n a . a.\ 

_ I ) 1 + 2tt i (1 -ji) 

If the angle of attack (a + 0 ) , referred to the axis of 

< 

zero lift, is expressed in degrees, we have 
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C = 8tt — 1ft + ' P)° (171 ) 

I 57.3 ^ - ; 

and thus 

C Z(A) = A^ A ) (a + B) C 



with 

»o _ ott a 1 tt 
A (.A) " 8 I 5773 4 



1 + f 1 



(27») 



1 + StrJ (l _|^) 
(A) 



The expression for the moment coefficient O mf A x remains 
unchanged and we have 

iO 



°»U) 8 *F 57.3 + T ° z Ur (2 ° } 



With the simplifying hypothesis of an elliptical lift dis- 
tribution, which assumes a x = 0, we would have 

C Z /a^ = O z , (37'") 

2 (A)E z (») 4 x + 2TT a v - ' 

L 

Rem ark I ,- The formulas relating to the expression of the 
mean lift coefficient Oz(^) for a win S of finite span depend 
essentially on the accepted law of lift distribution along the 
span. 

In using the approximate distribution represented by formula 
(26), we desired principally to attract the attention of tech- 
nical people to the necessity of considering, for rectangular 
wings, a distribution derivable from the aspect ratio. Conse- 
quently, we can employ the theoretical solution: established by 
Betz. 
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The formulas relating to the expressions for the moment 
coefficients 0 m(A) and C mQ , as well as the determination-, of 
the angle of zero lift, do not depend on the above law. of dis- 
tribution,, nor on the theory of finite span. These aerodynamic 
characteristics are inherent in the profile form, according to 
the theory of infinite span. 

Remark II.- The experiments conducted in the Aerodynamic 
Laboratory show that the mean lift coefficients experimentally 
determined are generally a little smaller . than, the theoretical, 
values computed by the preceding. formulas. This slight discrep- 
ancy results: first, from the approximation made in regard to 
the lift distribution; secondly, from the divergence between 
the real and the theoretical air flow; thirdly, from the ap- 
proximations inherent in the theory of finite span. 

Likewise, the experimental values of C mo and of £ are a 
little smaller than their theoretical values.- These differences 
seem to be due almost exclusively to divergencies between the 
real and theoretical flows. 

The comparison of the experimental and theoretical values 
of °z' c m 0 > and. (3 will permit figuring the discrepancies be- 
tween the experimental and theoretical aerodynamic characteris- 
tics for different groups of profiles. 

The induced polar is computed with a sufficient practical 
approximation by the formulas of Prandtl's theory. 
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2. Case of wings of elliptical or similar plan form. In 
this case, we still have a 1 = 0, but the expression for the 
mean lift coefficient becomes 



= C, 



or again 



Z (A)E ~ 2 (») ! + 21T ao 



° 2 (A)E = A ^ A)E ((X 



with 



l (A)E 



(38') 



= 8tt 



a r 



* 0 57 ' 3 1 + 3TT^> 

h 



a 0 and l 0 relating to the middle section. 

The expression for the moment coefficient C : 
calculated according to the plan form 



m (A) 



can oe 



Remarks.- The reduction factor to account for the discrep- 
ancy between theoretical and experimental lifts is a little 
smaller than the one relating to rectangular wings. 

3. «Case of wing oells .- This problem is more complex but, 
for a first approximation, we can use the preceding formulas by- 
giving to A the- value which is employed in calculations 

S 

of induced drag. 
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Diagrams of Aerodynamic Pressures on Wings of Finite Span 

The pressure coefficient at a point on any section of a 
finite wing is computed as already mentioned, taking care, how- 
ever, of the induced angle in the section under consideration. 
In other words, the direction! of the aerodynamic velocity, pass- 
ing through B 1 , makes with Ox an effective angle (a ~<P y ), 
cpy being the induced angle relative to the section located at 
a distance y from the middle section* 

For rectangular wings, having a distribution derived in 
terms of A, we have 

1. At the middle section: (y = 0) 

3TT » (1 _ Si) 
1 + 8n r (1 -jl) 

3. In any section located at the distance y from the 
middle section 




= 2tt St i _ |x + 3a 



x (i - »o) (30) 



or again 




1 + 2* 



in which i = a + 3 is the actual angle of attack with refer- 
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ence to the axis of zero lift. 

For elliptical or similar wings, with elliptical lift dis- 
tribution, we have a 0 

* - ■ rrS 1 (31) 

Aerodynamic Characteristics of Profile (Fig. 8) 

The chord length is laid out on the drawing, namely, 

I = 23.8 cm, 

while the approximate computation gives 

I = .2 (\ + OC ) = 23.54 cm 

The focus F falls on MF, making an angle 2 (7 - 8) = 
-4° with the first axis; and we have 

MF = = |%t = 5.38 cm 

We then measure 

AF = 6.15 cm, 
from which 

f-!£f -»•*»• 

The aerodynamic characteristics of the profile will then 

be 

0 Z/ = 8tt | (a + P ) = 6.8 (a + B) 
or, for (a + B) expressed in degrees, 
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V) ■ Ms (a + p)0 ■ °- 119 (a + p)0 - 

C mo = 0m (F) = 8n^i| ? = ? Vl° = 0.0536 

and ° m (A) = °^ 0536 + °» 258 c z(oo) # 

Lastly, for the experimental comparisons, we take the angle 
a c = - 1.3° between the tangent chord line and the axis Ox, 
From this it follows that the theoretical angle of zero lift, 
with respect to this chord, ha;s the value 

5.3° - 1.3° = 4° . 

Experimental Verification of the Aerodynamic Characteristics 

With the preceding profile drawing, we have dealt with ' 
rectangular wiigs with spans corresponding to the aspect ratios 
2, 3, 4, 5, and 6, with the same chord length 

I = 23,8 cm. 

The graph (Fig* 16) represents the experimental results 
of tests with these wings. 

We see that the experimental angle of zero lift is the 
same for all the aspect ratios and that its value is 3.6°, while 
the theoretical value was found to be '4°. 

The experimental C mQ , which is also the same for all as- 
pect ratios, is 0.045, while the theoretical value is 0.0556. 
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The mean slope of the curve (C Z) C m ) is 0.255 in place of the 
theoretical value of 0.258. 

Figure 17 represents the lift coefficients for the various 
aspect ratios converted to the infinite aspect ratio by the ap- 
plication of formula (27). The representative points are grouped 
along a mean line near the straight line representing the theo- 
retical lift of the profile. 

C z/ . = 0.119 (a + p).. 

The reduction coefficient k for the experimental and theo- 
retical lifts and angles £ is thus found to be k = 0.94. 

Lastly, Figure 18 represents the experimental points of the 
various polars converted to the single aspect ratio CA = 5. It 
is seen that the points thus corrected are grouped in a single 
polar blending with the experimental polar for A = 5. 

Experimental Verification of the Pressures 
on the Middle Section 

The pressure diagrams for the middle section were calcu- 
lated as already mentioned. The comparison of these theoretical 
pressures with the experimental pressures shows a very good agree- 
ment for the angles of* at tack at which the flow does not envelop 
important separations. 

This comparison between the theoretical and experimental 
aerodynamic characteristics confers an undeniable scientific. 
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character upon the investigation of wings with theoretical pro- 
files. Under these conditions, the experimental data, which, 
could he established by systematic tests of different profile 
groups,,, would be susceptible of a general interpretation* 

For example, in determining the law of evolution for the 
reduction coefficient between the theoretical and experimental 
characteristics for a group of theoretical profiles, we could, 
by a legitimate interpolation:, directly apply the suitable coef- 
ficient to all profiles of the group, without having to resort 
to a new experimental investigation. 

These remarks also apply to the aerodynamic characteristics 
which the present theories do not enable us to calculate. In 
particular, we could establish, for each group of theoretical 
profiles, the law of evolution: of the maximum C 2 and of the 
minimum- C x in terms of one or several characteristics of the 
drawing*' 

Thus, the choice of the profile to give an airplane wing, 
will be found to be completely determined by satisfying the re- 
quired conditions for the aerodynamic characteristics and struc- 
tural considerations. 

Translation by 

National Advisory Committee 

for Aeronautics. 
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Fig. 15 
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Figs j 12; 13, 14 
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Fig. 12 




Fig,i4 
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Fig. 17 
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Fig. 18 Lift coefficient converted to infinite aspect ratio. 

S.C. 281 (Von Mises) airfoil. Aspect ratios 6,5,4,3&2. 



